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1. Introduction 

The main goal of ultra-relativistic lieavy ion physics is to test the properties of matter, 
produced in nucleus-nucleus collisions at the highest energy densities accessible in the 
laboratory, and to develop an understanding of these properties from first principles of 
the fundamental theory of strong interactions. Quantum Chromodynamics (QCD). To 
this end, particle production in heavy ion collisions and its dependence on the formation 
of a dense system can be studied experimentally as a function of a large number of 
variables. These variables include the kinematic ones (such as centre-of-mass energy, 
transverse momentum and rapidity), as well as variables specific to heavy ion collisions, 
which typically control the size and shape of the collision region (such as the impact 
parameter or the nuclear number A of the colliding nuclei). 

Over the last two decades, relativistic heavy ion collisions have been studied 
experimentally at increasingly higher centre-of-mass energies at the Brookhaven 
Alternating Gradient Synchrotron AGS (y/s^ < 5 GeV), the CERN Super Proton 
Synchrotron SPS {y/s^ < 20 GeV) and the Brookhaven Relativistic Heavy Ion Collider 
RHIC {y/s^j^ < 200 GeV). As discussed in this article, the data collected in these 
experimental campaigns display remarkable generic trends as a function of system size 
and kinematic variables. 

The Large Hadron Collider LHC at CERN will study heavy ion collisions at a 
centre-of-mass energy ^/s^ = 5.5 TeV, which is almost a factor 30 higher than the 
maximal collision energy at RHIC. There has been a lot of work in recent years on 
benchmarking model calculations of heavy ion collisions to RHIC data and extrapolating 
them to the higher LHC energie^. The present topical review does not aim at such a 
comprehensive summary. Rather, the aim of this review is to identify the generic trends 
in the existing data and to discuss the consequences for our understanding of heavy ion 
collisions if these trends should persist or should fail to persist at the LHC. For those 
classes of measurements, where guidance from existing data is scarce (for instance for 
measurements at high-p^ or forward rapidity, where LHC is unique), we shall focus on 
generic features in the current model calculations, and discuss how they are expected 
to manifest themselves at the LHC. 

We believe that the identification of generic trends in the data and their agnostic 
extrapolation to LHC energies may help to sharpen our view on what is expected at 
LHC energies and what constitutes a surprise. In addition, many of the generic trends 
listed in this article did not yet find a fully satisfactory explanation. If they persist 
at the LHC, this would indicate that they should not be discarded as mere numerical 
coincidences, but should find an explanation in a future, more complete theory. 

I A comprehensive update on these efforts wiU be given in the proceedings of a recent CERN Theory 
Institute workshop [1]. 
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2. Multiplicity distributions 

In e~^e^ ^ gg — > X, where the partons produced initially carry perturbatively 
high virtuality, main characteristics of the longitudinal and transverse multiplicity 
distributions can be understood quantitatively from the dynamics of the perturbative 
parton shower (see e.g. [2]), despite uncertainties in the modelling of hadronization. 
In contrast, in hadronic collisions, event multiplicities and multiplicity distributions 
are dominated by processes involving non-perturbatively small momentum transfers; 
there are many models but an understanding of multiplicity distributions based on first 
principles is missing. Even in proton-proton collisions, the extrapolation of the charged 
particle multiplicity per unit rapidity dA^ch/dy from the Tevatron {^/s^ = 1.8 TeV) 
to the LHC (i/sj^ = 14 TeV) leads to results which vary by a factor 2 for models 
successful up to Tevatron energies [3]. For nucleus- nucleus collisions, the uncertainties 
in the predictions of minimum bias event multiplicities are of comparable magnitude. 
For instance, prior to the start-up of RHIC, model extrapolations of dNch/dy from the 
SPS centre-of-mass energy for Pb-Pb collisions {^/s^ = 17 GeV) to Au-Au collisions 
at RHIC {y/s^^ = 200 GeV) varied by a factor 2 approximately [4,5]. RHIC data lie at 
the lower end of the predicted range. This in turn has narrowed the range of predictions 
for the LHC. 

The lack of a fundamental understanding of multiparticle production in hadronic 
collisions is in marked contrast to several characteristic features, which persist over many 
orders of magnitude in y/s^ [6]: 

(i) Extended longitudinal scaling (limiting fragmentation [7]). 

Pseudorapidity distribution4§|, plotted in the rest frame of one of the colliding 
hadrons, fall on a universal, energy- independent limiting curve in the projectile 
fragmentation region. The region within which this limiting curve is valid, increases 
with energy, see figured! This is in contrast to the expectation that at high energies 
a boost-invariant plateau would develop around mid-rapidity. 

(ii) Factorization of ^/s^ and centrality/A-dependence. 

For all processes at a given centre-of-mass energy, the pseudorapidity distribution 
is the same basic distribution adjusted for the number of participants in the two 
colliding systems. The ^/s^- and A^part-dependences of dN^^/drj factorize. 

The apparent universality of these observations motivates an agnostic extrapolation 
to the LHC: If A-A data at LHC follow the same limiting fragmentation curve and if 
the trapezoidal shape of pseudo-rapidity distributions persists, then one expects for 
dA^^\/d?7 at LHC the solid line in figure [H This curve implies dN^p^/drj ~ 1100. 

§ Pseudorapidity 77 = tanhT^ pi/p is for many purposes a good approximation of rapidity y = 
ta,nh~^ pi/E. However, since multiplicity distributions are dominated by particles at small transverse 
momentum, there are visible differences: dN^^'^/drj is of trapezoidal shape (see figure [1]), while 
dN^/dy is of Gaussian shape. So, figure [T] does not imply a rapidity plateau. For A-A collisions, 
multiplying dN/drj at 77 = by the conversion factor « l.f turns out to be a good estimate for dN/dy 
at 2/ = 0. 
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Figure 1. Pseudorapidity distribution of cliarged particle production in Au-Au 
collisions at different centre-of-mass energies. Data are plotted in the rest frame 
of one of the colliding nuclei (full symbols), and mirrored at LHC mid-rapidity 
(open symbols). Agnostic extrapolations to the LHC are based on assuming limiting 
fragmentation and i) the saturation ansatz ([1]) (dashed line), or ii) a self-similar 
trapezoidal shape of the multiplicity distribution (solid line). Data from [8,9]. 



More generally, the persistence of extended longitudinal scaling in the high energy limit 
implies that at mid-rapidity, dN / dy can grow at most logarithmically with ^/s, except 
if there is a novel mechanism due to which the curvature of dN/ dy changes its sign twice 
between y = and the fragmentation region. 

Essentially all models of multiplicity distributions predict a power-law increase 
with y^s^. This is a rather generic consequence of perturbative particle production 
mechanisms, which become more important with increasing ^/s. However, the power- 
law dependence of naive perturbative implementations is too strong to be reconciled 
with RHIC data — this is arguably the main lesson learnt from the failure of many 
models at RHIC. Saturation models have received much attention recently, since they 
offer a fundamental reason for the very weak -^i- dependence of event multiplicities, 
namely the taming of the perturbative rise due to density-dependent non-linear parton 
evolution. Still, saturation models assume that multiplicity distributions at ultra- 
relativistic energies are calculable within perturbation theory, since they are governed by 
a perturbatively high, a/s- and A-dependent momentum (saturation) scale Q^^^ ^ oc . 
They predict essentially, that multiplicities at mid-rapidity rise oc Q^^^ ^ times transverse 
area. This leads e.g. to the pocket formula [10] 



2 diV^t^ 



= No ^/I^^^WGeV] N^ltt • (1) 



-/Vpart dr] 

Here A'^o = 0.47 is fixed by fitting to RHIC multiplicity distributions. The factors 
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A = 0.288 and 5 = 0.8 are constrained by fitting data on eA inelastic scattering. 
In this sense, the ^/s- and ^-dependences of ([1]) are predicted. The pocket formula 
([T]) provides an exphcit reahzation of the factorization property (ii) stated above, 
and accounts satisfactorily for the a/s- and A-dependences of charged multiplicity 
distributions from SPS to maximal RHIC energies [10]. For central Pb-Pb collisions at 
the LHC, equation ([1]) leads to dN^p^/drj ~ 1650, which corresponds to the maximum 
of the dashed line in figure [H An alternative model implementation of saturation physics 
ideas arrives at dN^p^/drj ~ 2200 [11]. A similar value is also predicted by the EKRT 
final state saturation model [12], which arrives at a ys^^Tj^'^^-dependence of the charged 
multiplicity at mid-rapidity. To arrive at a In y^s^- dependence in such schemes, one 
would have to invoke a mechanism, which amputates brutally the power-law tail in the 
spectrum for > (see e.g. the discussion of D. Kharzeev in reference [1]). 

A complete list of model predictions is beyond the scope of this article. We 
emphasize, however, that none of these predictions can be reconciled with the 
assumption that the so-far universal extended longitudinal scaling persists at LHC 
energies. As seen from figured], a distribution matching dN^py^/drj ~ 2200 at mid- 
rapidity has to fall off significantly steeper than what is consistent with limiting 
fragmentation. Also, the prediction dN^py^/drj ~ 1650, while being close to the 
maximum of what is consistent with limiting fragmentation, appears to deviate 
characteristically from the trapezoidal shape of all pseudo-rapidity distributions of 
charged multiplicity measured so far. 

To summarize and generalize this discussion: Either, the apparently universal 
"structure" seen in multiparticle production data at lower energies [6] is violated at 
the LHC. Then, this violation is likely to provide highly discriminatory constraints on 
the dynamics underlying multiparticle production. Or, the naive extrapolations of this 
structure to the LHC are confirmed. Then, the central djTiamical ideas advocated as 
explanations for the tamed growth of multiplicities up to RHIC energies will have to be 
revisited. So, starting with the first day of operation, data from the LHC are likely to 
have profound consequences for our understanding of the matter produced in nucleus- 
nucleus collisions at the LHC and at RHIC. 

3. Hadrochemistry 

The relative abundance of identified hadron species in heavy ion collisions follows a 
statistical ( "thermal" ) distribution pattern over a very broad energy range from SIS/GSI 
(\/^ ~ 2 GeV) up to RHIC (^/s^ = 200 GeV) [13]. For particle species, for 
which global constraints (such as total charge or flavour conservation) are statistically 
unimportant because of sufficiently large event multiplicities, particle ratios are well 
described by the grand canonical ensemble of a hadron resonance gas. The only 
parameters of this model are the temperature T at chemical decoupling, and the baryon 
chemical potential fiB- Fits to particle ratios reveal a characteristic ^s^-dependence 
of these parameters, see figure [2J The baryon chemical potential decreases by almost an 
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Figure 2. Thermal model fits at mid-rapidity of the hadrochemical freeze-out 
temperature T and the baryochemical potential as a function of the centre-of- 
mass energy ^s„„. Data points up to RHIC energies are taken from reference [14]. 
Data points at ^/s^ = 5.5 TeV are based on simple extrapolations of the observed 
tendency. 



order of magnitude from /i^ ~ 250 MeV at the SPS to jjB ~ 20-40 MeV at RHIC 
and is generally assumed to be very small {fiB ^ 10 MeV) at LHC mid-rapidity. 
This reflects the expectation that due to the large difference between projectile and 
mid-rapidity at the LHC, the mid-rapidity region is almost net-baryon free. The 
chemical decoupling temperature approaches a y^sj^-independent limiting temperature 
T ~ 160 — 170 MeV, which seems to be almost reached at RHIC and which is expected 
to persist at the LHC. These agnostic extrapolations to LHC energies have been used 
to predict the ratios of more than twenty hadron species, and there are computer codes 
implementing these model assumptions [15-18]. Early works on hadronic abundances 
used a strangeness saturation parameter [19] to account for enhanced strangeness 
production with increasing y/s^. This followed the idea that with increasing y/s^ 
a gluon-rich initial system is formed, which is more efficient in producing strangeness. 
In contrast, the above-mentioned models do not fix the strangeness content with an 
additional parameter. These models explain strangeness increase by a suppression of 
strange particles in low-multiplicity systems created at lower y'sj^. This suppression 
is due to exact conservation laws, whose implementation leads to deviations from 
the grand-canonical limit [13]. The resulting canonical suppression factors are known 
analytically. 

At RHIC and at the LHC, canonical suppression is unimportant for strangeness, 
which is produced abundantly. On the other hand, canonical suppression is important 
for open charm [20] and bottom production at the LHC. In this general sense, charm 
and bottom are likely to play at collider energies a similar role as strangeness played 
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during the fixed-target era of lieavy-ion pliysics. Tliere is, liowever, one important 
difference: because of tlieir large mass and small production cross section at tfiermal 
energies, tlie tliermal production of cliarm and bottom quarks is disfavoured. Tlius, 
except for the very first proposals of thermal charm production [21], models assume 
normally a production of heavy-flavoured quarks at perturbative rates, followed by a 
hadrochemical distribution of heavy-flavoured resonances according to the statistical 
model of hadroproduction [20] . In this approach, the centrality dependence of charm- 
fiavoured particle ratios tests the transition from a canonical to a grand-canonical 
description. 

The models discussed so far implement the idea of statistical hadronization, but 
they do not specify the dynamics leading to the hadronic final state. In particular, 
although the hadrochemical freeze-out temperature and baryochemical potential shown 
in figure [2] match within errors the QCD phase-space boundary determined in 
lattice QCD, the question whether this agreement implies the existence of dynamical 
thermalization processes lacks a more detailed support [22] . We conclude this section by 
listing three possibilities of how heavy- ion collisions at the LHC may help to elucidate the 
microscopic dynamics underlying hadronic abundances: First, it has been conjectured 
that charmonium production at the LHC is strongly enhanced above the perturbatively 
expected rates due to the recombination of c and c quarks originating in distinct hard 
partonic interactions [23] . Establishing such a novel charmonium production mechanism 
could provide a strong indication that thermalization processes affect hadrochemical 
distributions. This is so, since transport of exogamous c- and c-partners is likely 
to be necessary to make their coalescence possible. Second, we mention in passing 
a proposal that the abundances of strange hadrons could exceed significantly grand- 
canonical predictions. This is seen in models in which a sudden hadronization of an 
equilibrated plasma leads to strangeness over-saturation in the hadronic phase [24]. 
While arguably speculative, this example illustrates that the extrapolation of (standard) 
statistical model predictions to the LHC can serve as a powerful baseline on top of 
which novel dynamical effects may be established. Third, the masses and widths of 
the resonances entering the statistical operator of a resonance gas model may receive 
significant medium-modifications. This may lead to deviations of some broad resonances 
(such as p-mesons) from the grand canonical ensemble. While not specific for the LHC, it 
is likely that any confirmation and extended systematics of such deviations will provide 
constraints on the in- medium dynamics of hadronic resonances near freeze-out. 

4. Transverse-momentum spectra at low 

The transverse momentum dependence of identified single inclusive hadron spectra has 
been studied in heavy ion collisions from SIS/GSI up to RHIC energies. Here, we follow 
the common practice of discussing these spectra as a function of their transverse mass 
niT = \/m? -(- p\^ rather than their transverse momentum. The m^-distributions are 
commonly characterized by a mono-exponential fit dA^/dm^ oc m^''^ exp [— m^/Tinv] [25], 



Predictions for the LHC heavy ion programme 



8 



where the inverse-slope parameter Tiny characterizes the steepness of the distributions. 
The fit value Ti^y can depend on the m-r-range over which the fit is performed, but 
despite the ensuing uncertainties, the following generic features can be stated: 

(i) Linear mass-dependence of Tjnv 

At given ^/s^, the inverse-slope parameters of pions, kaons, and protons increase 
approximately linearly with the particle rest mass. In particular Ti^^^ ^± < Ti^^^ ^± < 
Tinv.p- The hadrons, which do not follow this systematics, such a the multistrange 
S, the Q or the J/^', show roughly the same value Tjnv A seemingly common 
denominator of these latter hadrons is their relatively small interaction cross 
sections with the expected constituents of the medium [26,27]. 

(ii) Increase of Tlnv with ^/s^. 

The spectra become flatter with increasing ^/s^, so the inverse slope parameter 
Tinv increases. Also, the mass-dependence of Tlnv for pions, kaons and protons 
increases with ^/s^ [26]. 

The inverse-slope parameter Tiny has been interpreted as a blue-shifted temperature, 
resulting from a combination of thermal emission from a source of temperature Tf.o. 
at freeze-out, and the collective motion with average transverse velocity (Pt) of this 
particle emitting source. While Tf o and (/3t) are difficult to disentangle on the basis 
of single inclusive spectra alone, two-particle correlations may be used to separate both 
contributions [28]. The increase of T^^^ with particle rest mass [point (i)] is consistent 
with the interpretation of Tiny as a blue-shifted temperature [25,29], if the particles 
decouple at the time of freeze-out. On the other hand, particles which decouple earher 
due to their smaller cross sections would acquire less radial flow, and thus have a smaller 
Tiny. In this radial-flow picture, the growing difference between inverse-slope parameters 
as ^/s^ increases [point (ii)] is ascribed to the growth in the average transverse velocity 
{/3t)- 

The above features emerge naturally in fluid-dynamics models without [30, 31] or 
with [32] viscous corrections. We are not aware of full fluid-dynamic studies of these 
spectra at LHC energy. However, there arc extrapolations of the fit parameters Tf.o. and 
{Pt), extracted from the approximately exponential mr-spectra [26]. In these studies, 
the kinetic freeze-out temperature saturates at a value Tf.o. ~ 120 MeV at SPS energy, 
which may be expected to persist up to the LHC. The parameter {Pt) is found to 
increase with ^/s^^, it reaches ~ 0.55 at RHIC, but it remains unclear whether 
this increase is smooth [26] . 

The statements made here are obtained within a blast- wave model [26,29] or within 
a fiuid-dynamic picture supporting such a model. They generally indicate an increase of 
Tinv with y/s^. Here, the real physical issue is to assess whether it is really a collective 
hydrodynamic mechanism, which determines the energy-dependence of Tiny. Since the 
slope of spectra evolves significantly with ^/s^ also for more elementary proton-proton 
collisions, in which flow is either absent or at least different, it is not straightforward to 
disentangle flow effects from other mechanisms. This makes quantitative comparisons 
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difficult and leads to a lack of precision in extrapolations to the LHC. We now turn to 
the azimuthal dependence of single inclusive transverse momentum spectra, where these 
difficulties appear to be less severe. 



5. Azimuthal anisotropy in \aw-px particle production 

In non-central collisions, the impact parameter selects a preferred direction in 
the transverse plane, breaking the rotational symmetry around the beam axis. 
Consequently, observables can depend on azimuth, measured with respect to the reaction 
plane spanned by the impact parameter and the beam direction. The azimuthally- 
dependent effect most widely studied both experimentally, over the whole range of 
available collision energies, and theoretically, is the azimuthal anisotropy in particle 
production, often referred to as "anisotropic (transverse) flow" [31]. A non- vanishing 
anisotropic flow exists only if the particles measured in the final state depend not only 
on the physical conditions realized locally at their production point, but if particle 
production does also depend on the global event geometry. In a relativistic local 
theory, this non-local information can only emerge as a collective effect, requiring 
interactions between many degrees of freedom, localized at different points in the 
collision region. In this general sense, anisotropic flow is a particularly unambiguous 
and strong manifestation of collective dynamics in heavy-ion collisions. Here, we discuss 
its features for hadrons at low-pr- 

Azimuthal anisotropics in particle production are most conveniently characterized 
by performing a Fourier expansion of the single-particle distribution 

-IT^ = TT A [1 + - + cos 2(0 - + ■ ■ ■] , (2) 

d^p^ ay 271 pt dpr dy 

where denotes the azimuth (in the laboratory frame) of the reaction plane. The 
coefficients f„ = (cosn(0 — $i?)), where angular brackets denote an average over many 
particles and events, quantify the asymmetry. These coefficients are studied at all 
available energies as a function of transverse momentum, (pseudo)rapidity, centrality of 
the collision and for various identified particle species. 

"Elliptic flow" , the second Fourier coefficient V2 is the best studied one. A positive 
(resp. negative) value of f 2 indicates an excess of particle production in (resp. orthogonal 
to) the reaction plane. The dependence of V2 on centre-of-mass energy is known over 
three orders of magnitude, see figure [31 It can be understood qualitatively in terms of 
the following simple picture of a collective dynamics: At ^/s^ < 2 GeV, the incoming 
nuclei transfer angular momentum to the nuclear matter in the overlap zone. The fast- 
rotating "compound nucleus" thus formed emits fragments, which lie preferentially in 
the reaction plane {v2 > 0). As the centre-of-mass energy increases, the nuclear matter 
in the almond-shaped overlap region of the incoming nuclei is increasingly compressed 
in the collision. However, the parts of the nuclei that lie outside this overlap region 
("spectators") block the way for the compressed matter to expand within the reaction 
plane. They squeeze-out the compressed matter orthogonal to the reaction plane 
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Figure 3. y^s^^-excitation function of i'2(2/ = 0) in mid-central collisions. Data are 
taken from the compilation in reference [33]. 



(f2 < 0). Further increasing y/s]^, the spectators are then fast enough to free the 
way, leaving behind at mid-rapidity an almond-shaped azimuthally asymmetric region 
of dense QCD matter. This spatial asymmetry implies unequal pressure gradients in 
the transverse plane, with a larger gradient in the reaction plane ("in-plane") than 
perpendicular to it. As a consequence of the subsequent multiple interaction between 
many degrees of freedom, this spatial asymmetry leads to an anisotropy in momentum 
space: the final particle transverse momenta are more likely to be in-plane than "out- 
of-plane", hence f2 > 0, as predicted in [34]. 

The momentum space asymmetries measured at collider energies are relatively 
large. Since the prefactor of the cosine term in equation ([2]) is 2v2, a pr-averaged value 
V2 = 0.05 corresponds to a 20% variation of the average particle yield as a function of 
the angle with respect to the reaction plane. At high p^, where second harmonics at 
RHIC approached values as large as V2 = 0.2, there are more than twice the number of 
particles emitted in the reaction plane than out-of-plane. Elliptic flow is an abundant 
and very strong manifestation of collectivity, which shows remarkable generic trends: 

(i) The pT-integrated f2(?7) shows extended longitudinal scaling [35]. 

In contrast to dA^/d?7, V2{ri) is not trapezoidal but triangular, see figure As 
seen clearly from figure HI longitudinal scaling of pr-integrated f 2 persists up to 
mid-rapidity. 

(ii) The pT-shape of the charged-hadron V2 has a characteristic breaking point. 

At transverse momenta below pt — 2 GeV/c, where data are known from SPS 
and RHIC, V2 is found to have an approximately linear rise with pt- Around 
Pt — 2 GeV/c, this rise levels off rather abruptly. The energy-dependence of this 

II The pT-averaged value of V2 is dominated by values of the transverse momentum close to (pt), so 
that V2{ri) and V2{y) are similar, in contrast to dA^/d?7 and dN/dy. 
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Figure 4. The elliptic flow V2, averaged over centrality (0%-40%), at various collision 
energies. Data (full symbols) from PHOBOS [35] and STAR [36] are plotted as a 
function of 77 — j/beam and reflected (open symbols) across the LHC — j/boam value. 



pr-shape is not fully clarified: At low pt, the slope of V2{pt) was reported to rise 
either slightly [37] or significantly [38] across SPS and RHIC energies. Also, it was 
reported [38] that the slope of V2{pt) saturates at RHIC energies and is essentially 
constant between ^/s^ = 62.4 GeV and 200 GeV. In this case, the increase of 
the pr-averaged f 2 would be entirely due to the increase in the mean transverse 
momentum of particles. 

(iii) The p^-dependent V2 of identified hadrons shows mass-ordering at small pr and 
displays a constituent-quark counting rule at intermediate px- 
For a fixed, sufficiently low transverse momentum, SPS and RHIC data show 
generically that V2{pt) decreases with increasing mass of the particle species. Above 
a critical pt ~ 1-5 GeV/c, mass-ordering ceases to be valid and V2{pt) follows to 
a good approximation a simple quark counting rule, namely that V2{pT l^q) /uq is 
roughly independent of the particle species [39]. 

What are the implications if these trends persist or do not persist at the LHC? 
First, if longitudinal scaling of V2 persists, then f 2(?7) grows proportional to In y/s^. In 
this case, one expects ^2(^7 = 0) ~ 0.075 for Pb-Pb collisions in mid-central collisions. 
This follows from the extrapolations, shown in figures [3] and |H To the best of our 
knowledge, neither the triangular shape of the rapidity dependence of V2, nor the 
approximately linear In y^s^-dependence emerges as a natural consequence of existing 
dynamical models. In particular, extrapolating models of ideal hydrodynamics from 
RHIC to the LHC, one arrives at values not exceeding V2{ri = 0) ^ 0.06 for event 
multiplicities shown in figure [1] [30]. Also, the proportionality V2{y) oc dN/dy does not 
hold in models presupposing local equilibrium (i.e. the formation of an almost perfect 
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fluid). If tlie matter produced at RHIC mid-rapidity is in local equilibrium, then one 
expects at the LHC deviations from the triangular shape of ^2(^7) in an extended region 
~3 ^ y ^ 3, around mid-rapidit}{^. In summary: extrapolations to the LHC of the main 
dynamical mechanism advocated to underlie elliptic flow at RHIC (namely perfect fluid 
dynamics) are at odds with extrapolations to the LHC of the generic trends observed in 
elliptic flow measurements up to RHIC energy. As a consequence, establishing whether 
these trends persist at the LHC provides a novel independent test for our understanding 
of the properties of matter at the LHC and at RHIC. 

Second, comparing measurements of V2{pt) from LHC and RHIC at low px will 
finally allow us to establish to what extent the pr-slope of V2{pt) changes with y/s^. 
This is of interest, since a ysj^-independent slope of V2{pt) would imply e.g. that 
the increase of the pt- integrated V2{ri) with y/s^ arises solely from the increase of 
the average transverse momentum. Existing dynamical models of V2 do not invoke 
the increase of the rms transverse momentum \/Jj^) with ^/s^ observed in hadronic 
collisions; so, establishing a major role of a/ {p^) in V2 may prompt significant revision 
in our interpretation of elliptic flow. Moreover, LHC data for V2{pt) at intermediate 
Pt will test to what extent the breaking point of V2{pt) depends on centre-of-mass 
energy, collision centrality, or choice of nuclei, while existing RHIC data hint at a 
very small sensitivity on these [42]. In the current discussion of RHIC data, one 
emphasizes that the approximately linear increase of V2{pt) at low pt is consistent 
with the transverse expansion of an almost perfect liquid [31], while the breaking point 
in V2{pt) arises from the rather abrupt onset of dissipative effects at higher pt [43]. 
However, transport models and dissipative hydrodynamics [43] can account for the linear 
increase of V2{pt) at small px as well, if initial conditions are chosen appropriately, and 
a detailed understanding of the dynamical origin of the breaking point is missing and 
may proflt from knowledge about its y^s^-dependence. 

Third, mass-ordered V2{pt) at mid-rapidity are predicted in transport ap- 
proaches [44,45] as well as in hydrodynamical models [30,31], but details of the predicted 
mass hierarchy vary between models. For instance, two different ideal-fluid dynamics 
approaches predict scaling laws with different variables: the decoupling of all particle 
species from the same collective-flow fleld is argued in reference [46] to lead to a common 
V2{pT/'m) for all hadrons; while in the Buda-Lund framework [47] a scaling of V2{pt) 
with the square of the transverse rapidity yr = |ln[(mr +pT)/{jn'r —pr)] is predicted. 
So, while mass ordering of V2{pt) for light hadrons is expected to persist qualitatively, 
its quantitative manifestation may help to differentiate between models of the collision 
dynamics. 

Moreover, at the LHC, the elliptic- flow parameters V2{pt) of D- and 5-mesons 
will provide yet another test of the mass-ordering of V2- Quantitative comparison to 
predictions of dynamical approaches — be it fluid dynamics, a Langevin description, or 
a transport model — have been argued to give insight on the possible thermalization [48] 

^ However, in models which presuppose incomplete local equilibration [40,41], one finds 1^2(2/) oc AN /dy, 
thus accounting for the triangular shape of V2{ri). 
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of heavy quarks, in particular their mean free path [49], and their hadronization 
mechanism [50,51]. Quite generally, heavy quarks in equilibrium have a larger V2{pt) 
than non-equilibrated ones; and the corresponding mesons have larger elliptic-flow values 
if they form through coalescence (involving a light quark) than if they come from heavy- 
quark fragmentation. 

Mass-ordering persists essentially up to the breaking point in the p^-shape. Above 
this point and up to ~ 5—6 GeV/c, a quark-counting rule was observed at RHIC for light 
hadrons, defining a region of intermediate transverse momenta. (New rig-scaling rules 
at low transverse momentum were also recently reported at RHIC [42,52], which to our 
knowledge have no theoretical explanation in the existing literature.) The persistence 
of such rules at LHC would constrain dynamical models, especially those covering both 
the low- and intermediate-pr regions. 

We now turn to the centrality dependence of V2. It has been suggested to classify 
finite impact-parameter collisions for nuclei of different size in terms of the surface 
5* of the transverse overlap region and the eccentricity e of this surface. (We note 
that fluctuations [53, 54] and uncertainties in the initial conditions [55, 56] make the 
specification of e and S somewhat model-dependent.) One observes that data of 
^2(1/ = 0)/e from AGS to RHIC energy show an apparently universal, linear dependence 
if plotted against {1 / S) dN^^ / dy (see e.g. figure 25 in reference [57] or figure 15 in 
reference [58]). Modelling the matter produced in heavy ion collisions in terms of a 
perfect fiuid, one finds that the above-mentioned linear dependence levels off above 
values of {1 / S) dN'^^ / dy corresponding to central RHIC energies. Thus, as in our 
discussion of several other classes of measurements, a naive extrapolation of these data 
on centrality from RHIC to LHC is at odds with a naive extension of the main dynamical 
explanation advocated to underlie elliptic flow V2- On the other hand, there are models 
which may account for a further increase of ^2/6 due to non-equilibrium phenomena in 
the initial state [41], or due to a change in the relative contributions of hadronic and 
partonic rescattering effects as a function of y/s^ [59], or due to a significant change 
in initial conditions [60]. 

The second most studied fiow harmonic is Vi, "directed fiow", which quantifies 
the average momentum acquired by the particles along the impact-parameter direction. 
Two generic trends can be identified in existing data: 

(i) The pT-integrated vi{y) is linear around mid-rapidity. 

At mid-rapidity, vi vanishes by symmetry in collisions between identical nuclei. A 
linear increase of Vi{y) around y = is observed across AGS and SPS energies. 
The slope dvi/dy decreases with increasing beam energy. 

(ii) The p^-integrated Vi{ri) shows extended longitudinal scaling [61]. 

Above SPS energy, one finds that Vi{ri) is positive in the "projectile" {rj > 0) 
fragmentation regioii^, then becomes negative for rj < yheam, reaches a minimum 

+ This is a choice, inherited from fixed-target studies at lower energies, where the bounce of the 
projectile off the target is taken to define the positive direction. The absolute sign of vi is not 
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for -2 < ?7 - ?/boam < -1, and increases. 

The requirement that Vi be zero at midrapidity imphes either a breakdown of the 
longitudinal-scahng property, or that Vi (rj) vanishes in an extended region around y = 
(so that the slope dvi/dy is also zero), the size of which increases with beam energy. 
In either case, the pr-integrated vi{y) at LHC will be smaller in absolute value than 
0.01 up to rapidities y ^ 4 — 5. A deviation from this expectation would indicate a 
physics effect that manifestly breaks Bjorken boost-invariance and which becomes more 
pronounced at higher y/s^. We are not aware of any suggestion of such an effect at 
these energies. Thus, testing a non-trivial dependence of vi at LHC and extending the 
RHIC systematics is likely to require measurements at far forward rapidity, which are 
experimentally challenging. 

We finally comment on the fourth anisotropic-flow harmonic V4. It has been argued 
that at given transverse momentum and rapidity, v^ipxiy) > \v2{pT,yY [41,46], the 
lower bound being attained if and only if the matter expands like a perfect fluid at the 
time when anisotropic flow develops. At RHIC, one has reported ^4/^2 — 1-2 [62]. If one 
takes this number as an indication of incomplete equilibration, and if one assumes that 
equilibration mechanisms are more efficient at the LHC, one predicts 0.5 < Vi/vl < 1.2 
at the LHC. Without such assumptions, one still observes that both RHIC data and 
model calculations lead us to expect a value of v^/vl of order unity. 

6. Femtoscopy 

The header femtoscopy summarizes a class of measurements that give access to the 
spatio-temporal extension and collective dynamics of the matter produced in heavy-ion 
collisions [63] . This includes in particular identical two-particle momentum correlations 
C(K, q) of relative pair momentum q and average pair momentum K. These are often 
analysed in the Bertsch-Pratt parametrization 



Here, the indices i,j label a Cartesian coordinate system with axes pointing along 
the longitudinal beam direction {longitudinal or "1"), parallel to the transverse pair 
momentum {out or "o") and the remaining {side or "s") one. The radius parameters 
i?jj(Kx, i^^i) combine information about the spatial and temporal extension of the 
particle-emitting source at freeze-out. They do not measure the extension of the entire 
collision region, but they measure the generally smaller "homogeneity regions" , i.e. the 
part of the source radiating particle pairs with pair momentum (K^, Kl). Data of these 
HBT radii, taken at the SPS and RHIC, display several generic trends: 

(i) Almost linear scaling with (dA^ch/d?])^^'^. 

As seen in figure El the diagonal radius parameters in the out-, side-, and 

measurable, only its changes in sign are. 




(3) 
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Figure 5. The diagonal HBT radius parameters R\, Rs and Ro at mid-rapidity 
plotted versus event multiplicity for different centre-of-mass energies, including an 
extrapolation to LHC. Data from the STAR collaboration taken from reference [64]. 
Except for the 200 GeV Au-Au results, all data are preliminary, but consistent with 
previously observed trends [63] . 



longitudinal directions scale approximately linearly with the third root of the 
charged particle multiplicity per unit rapidity. The entire ^s^^^^ -dependence 
appears to arise via the dependence on (diVch/dr])^/^. 

(ii) Universal mr-dependence. 

Values for the HBT radius parameters of different particle species fall on a universal 
curve if plotted versus transverse mass, = ^/m'^ + K^. In particular, the ansatz 

R^iKr) ^ c,-^ , (4) 

provides a fair description of the universal mT^-dependence. Here, the parameters 
ai and Cj take the same universal values for all hadron species. Fit values of are 
~ 0.5, albeit with large variations. 

(iii) Rl ^ Rl 

The three diagonal radius parameters differ in size, but this difference is small 
compared to their absolute value. In particular, i?^ ~ _Rg. 

The extrapolation of the linear-scaling trend (i) from RHIC to the LHC leads 
to values for the HBT radius parameters that exceed significantly those measured 
previously, see figure [51 For instance, one finds i?i(150 MeV < Kt < 250 MeV) ~ 8 fm 
for dA^ch/d?7 = 1100, and even larger values for larger event multiplicities. For central 
collisions at mid-rapidity, the off-diagonal radius parameters Rij vanish by symmetry, 
and the product V = RyRg Ro provides a working definition of the spatial volume of the 
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homogeneity region. The hnear increase of all three HBT radii with (dAi^ch/d?])^/^ is then 
consistent with the statement, that hadrons freeze-out from the collision system at an 
universal phase-space density [65]. Changes in hadrochemical composition, transverse 
flow or temperature gradients may lead to deviations from this freeze-out criterion 
at universal phase-space density [66]. In particular, an increase of the fraction of 
baryons {p + p + higher resonances) at fixed density leads to an increase of HBT 
radius parameters, since baryonic cross sections are larger than mesonic ones and 
thus delay freeze-out [67]. In contrast, larger flow or temperature gradients tend to 
decrease the HBT radius parameters, since they narrow the spatial extension within 
which identical particle pairs show significant quantum-mechanical interference. The 
linear extrapolation Ri oc (dA^ch/d//)^^'^ provides an agnostic baseline on top of which 
dynamical changes may be established. 

In models with flow-dominated freeze-out scenarios, the mr-dependence of the 
transverse HBT radii steepens as flow increases and/or temperature decreases. However, 
the numerical significance of this effect is model-dependent and other factors (such as 
effects from resonance decay contributions or the opacity of the produced matter) may 
play a role as well. Even if one assumes that transverse flow effects increase with 
y/s^ from RHIC to LHC, the question whether the mr-dependence of transverse HBT 
radii steepens, appears to be too subtle to have a robust, model-independent answer. 
Moreover, in model studies, transverse flow may manifest itself in a f^T-dependence of 
transverse HBT radii, which cannot be absorbed in an mT-dependence. We are not 
aware of sufficiently generic LHC predictions for this interesting class of measurements. 
As an aside, we note that for establishing the trend displayed in figure [5l data for different 
centralities and centre-of-mass energies should be compared at the same transverse 
momentum — otherwise, the strong mT-dependence (jl]) may mask the apparently 
universal dependence on {dNdJdrjY^^. 

Under relatively mild model assumptions, the difference B^{Kt) — RI{Kt) can be 
related to the lifetime of the particle-emitting source. The smallness of this quantity 
has been dubbed puzzling, since many models of the source dynamics predicted large 
lifetimes (for instance as a consequence of a first order phase transition or rapid 
crossover) [68]. Also, fluid dynamics strongly over-predicts — R^, if it is not 
supplemented by extended hadronic scattering mechanisms [69], or by other significant 
modifications of the final state [70]. In the energy range from RHIC to LHC, no 
mechanism is expected to set in, which could modify the value R^{Kt) — RI{Kt) 
significantly. 

The analysis of femtoscopic information at LHC is not limited to the three generic 
features listed above. It includes the analysis of HBT radii with respect to the reaction 
plane [28] , the physics hidden in the intercept parameter A of equation ([3]) , and a rapidly 
growing field of non-identical particle correlations. So far, however, data on these classes 
of observables are too scarce to serve as a robust baseline for extrapolations to the LHC, 
and we refer to a recent review article for further discussions [63]. 
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7. Single inclusive high-pt spectra in A-A and pA collisions 

The study of single inclusive hadron spectra at high transverse momentum has lead 
to some of the major discoveries of the RHIC heavy ion program [71]. At the SPS, 
kinematic constraints limit the analysis of transverse momentum spectra in practice to 
Pt < 3 — 4 GeV/c. In contrast, the RHIC program studied at higher centre-of-mass 
energy various hadron species in a range up to px < 10 — 20 GeV/c, where perturbative 
production mechanisms are known to account for the single inclusive spectra in hadronic 
collisions. At the LHC, the transverse phase space accessible for such measurements 
increases by another factor ~ 10. The limited kinematic reach of the SPS prompts us 
in the following sections to base extrapolations to the LHC on RHIC data only. 

The nuclear modification factor R\q characterizes how the production of a hadron 
/i in a nucleus-nucleus collisions A-B differs from its production in an equivalent number 
of proton-proton collisions. 



Here, (-/V^u) average number of inelastic nucleon-nucleon collisions in a given 

centrality class. This number is typically determined in a Glauber-type calculation. 
The nuclear modification factor depends in general on the transverse momentum pt and 
pseudo-rapidity rj of the particle, the particle identity h, the centrality of the collision 
and the orientation of the particle trajectory with respect to the reaction plane (which 
is often averaged over). In the absence of medium effects, — 1. 

7.1. The nuclear modification factor at mid-rapidity 
RHIC data on Raa show the following generic features: 

(i) Characteristic centrality dependence of Raa and i?dAu- 

For the most peripheral centrality bin, the nuclear modification factors measured 
at RHIC arc consistent with the absence of medium-effects in both nucleus-nucleus 
{Raa ~ 1) and deuterium-nucleus (i?dAu ~ 1) collisions [72-75]. With increasing 
centrality, Raa decreases monotonically. In d-Au collisions, the opposite centrality 
dependence is observed with maximal values RdAu ~ 1-5 around pt = 3 — 5 GeV/c 
in the most central bin. Accordingly, one observes that Raa depends on the 
azimuth with respect to the reaction plane, with a smaller Raa out-of-plane [76] — 
equivalently, V2(pt) is positive at high px [36]. 

(ii) Strong and apparently pT-indcpcndcnt suppression of Raa at high Pt ■ 

In y/s^ = 200 GeV, 5-10% central Au-Au collisions at mid-rapidity, one observes a 
suppression of high-p^ single inclusive hadron yields by a factor ~ 5, corresponding 
to R\uAu{Pt) — 0.2 for > 5 — 7 GeV/c. Within experimental errors, this 
suppression is pT-independent for higher transverse momenta in all centrality 
bins [77-79]. 



■AB 



{Pt, r], centrality) 




(5) 
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(iii) Independence of Raa on hadron identity. 

For transverse momenta Pt > 5 — 7 GeV/c, all identified (light-flavoured) hadron 
spectra show a quantitatively comparable degree of suppression. There is no 
particle-species dependence of the suppression pattern at high pt- 

(iv) The photon spectrum is consistent with perturbative expectations. 

For single inclusive photon spectra, the nuclear modification factor shows 
mild deviations from -RXuAu ~ -'- l^^]- Within errors, these are consistent 
with perturbative predictions taking into account the nuclear modifications of 
parton distribution functions (mainly the isospin difference between protons and 
nuclei) [81]. 

(v) Raa shows a characteristic baryon- meson difference at intermediate px- 

At intermediate pr, 3 GeV/c < pr < 6 GeV/c say, the nuclear modification 
factor for mesons is smaller than the one for baryons [82]. Within experimental 
uncertainties and irrespective of the hadron mass, all identified meson spectra show 
a similar degree of nuclear suppression, and so do all identified baryon spectra. 

For hadronic collisions, the perturbative QCD factorized formalism can account 
systematically for single inclusive hadron spectra at sufficiently high transverse 
momentum, by convoluting ("incoming") parton distribution functions, with hard, 
partonic scattering matrix elements, and with ("outgoing") parton fragmentation 
functions. In nucleus-nucleus collisions, one aims at identifying the leading medium- 
length enhanced nuclear effects which modify this factorized formafism [83]. In this 
context, the notions "ingoing" and "outgoing" become physically relevant, since the 
hard production process is placed within the spatio-temporal geometry of a nuclear 
collision. The zeroth order question is whether the dominant medium modification 
is accumulated during the incoming or outgoing stage of its prolonged interaction 
with the medium. Experimentally, this can be addressed by systematically varying 
the final state effects; for instance by varying the outgoing in-medium path length 
via centrality measurements, or by switching off final state effects by comparing A-A 
collisions with HA collisions. Theory addresses these dependencies in model studies, 
which supplement the perturbative QCD factorization approach to single inclusive 
hadron spectra with medium-modifications in the initial and final state, taking the 
spatio-temporal distribution of matter during the A-A coUision into account [84-86] . 

RHIC data prove that high-p^^ hadron suppression is predominantly a final state 
effect by establishing that the suppression is not seen in d-Au and that it increases in 
A-A with increasing centrality and thus with increasing in-medium path length in the 
final state [point (i)]. Moreover, the independence of R\j^ on hadron identity [point 
(ii)] at high pt gives support to the picture that the final state medium modification 
of parton fragmentation is of partonic nature, i.e. that it occurs prior to the onset of 
hadronization. This is so, since hadronic states would present absorption cross sections 
which can be expected to differ significantly with hadron identity, and should thus lead 
to a hadron-specific splitting of the nuclear suppression factor i?^^, which is not observed 
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above pr > 5 — 7 GeV/c. In addition, the single inclusive photon spectrum indicates that 
initial state effects at high pt are small and that they may be accounted for by nuclear 
modified parton distributions [point (iv)]. From these arguments, one concludes that 
the suppression of high-p^ single inclusive hadron spectra in nucleus-nucleus collisions 
is due to a partonic, medium-length dependent final state effect. 

From an agnostic point of view, one wonders whether the generic features (i)-(v) 
persist at the 30 times higher y/s^ explored at the LHC, and how these features evolve 
in the much wider transverse momentum range accessible at the LHC. To discuss these 
questions, we first turn in section [7. 1.11 to the extrapolation of models that describe the 
generic features observed in RHIC data. Then, we provide in section 17. 1.21 a list of effects 
that may become important at the LHC and could lead to characteristic deviations from 
the generic features observed at RHIC. 

7.1.1. Extrapolations of parton energy loss models to the LHC. The generic features 
(i)-(iv) are naturally accommodated in models which supplement the perturbative QCD 
factorization approach (using nuclear parton distribution functions) with a mechanism 
which degrades the energy of the leading outgoing partonic fragment due to its 
propagation in matter. Two classes of mechanisms have been explored in Feynman 
diagrammatic detail: coUisional [87] and radiative parton energy loss [88-92, 100]. 
Radiative energy loss, that is the medium-enhanced splitting of the energetic parton, 
is the dominant mechanism at high pt [84-86]. Up to which pt subleading collisional 
effects are numerically significant and how they could be disentangled from radiative 
ones is a matter of ongoing debate [93-95]. In the following, we limit our considerations 
to radiative energy loss, which gives a fair description of RHIC data above pr > 
7 GeV/c [96-99]. This may provide a baseline on top of which collisional contributions 
can be estabhshed. 

In radiative parton energy loss models, only one medium- dependent model 
parameter enters, the so-called jet quenching parameter g(r) (or a reparametrization 
of it), which depends on the time r after the collision. In model studies, it is often 
expressed in terms of the energy density e(r) [100] 

q{T)=ce'/\T), (6) 

where c is assumed to be a time- and temperature- independent constant. An estimate 
based on perturbatively weak interactions between the hard parton and the medium 
gives c ~ 2 [101]. In contrast, fitting parton energy loss models to RHIC data, several 
groups found much larger values, c > 8 [98,99,102]. Assuming a 1-dimensional Biorken 
expansion of the produced matter with g(r) = qoTo/r, this translates into a valufB 

g(r = 1 fm/c) > 4 GeVVfm . (7) 

* Since medium-induced gluon energy radiation in an expanding medium depends on the line-average 
q = dr r g(r) [103], several model studies quote q for an average in-medium path length [98,99]. 

However, q depends strongly on L, which differs for each parton. This introduces additional 
uncertainties. To bypass these problems, we quote here the value c, or equivalently qir ~ 1 fm/c) 
which can be unambiguously extracted from all studies quoting q. 
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The precise value of q consistent with RHIC data is currently debated, but it is 
generally thought that q is significantly larger than the perturbative estimate c = 2 
in reference [101]. 

The jet quenching parameter q has a rigorous field theoretical definition in terms 
of the short-distance behaviour of the target expectation value of a light-like Wilson 
loop [104]. For a class of non-Abelian thermal gauge field theories, which are known 
to have a gravity dual, non-perturbative evaluations of q have established recently that 
the seemingly small scale of an initial temperature of T{t = Ifm/c) ~ 300 MeV does 
indeed give rise to a jet quenching parameter g(r=lfm/c) numerically consistent with 
the apparently large lower bound of equation ([71) [104]. These studies also show that 
the ratio of jet quenching parameters q of different thermal field theories is determined 
by the square root of the ratio of their entropy densities [105]. In earlier studies, q was 
taken to be proportional to the event multiplicity, and large values up to (?lhc — '''^rhic 
were explored [98,99]. On the other hand, the multiplicity extrapolations to the LHC 
shown in figure [T] would indicate a much smaller value of ^lhc' — given that a 
multiplicity scaling of q is an additional model assumption — a very mild increase of 
say gLHC — l-25gRHic is conceivable. Since q is the only medium- sensitive parameter 
in a class of model studies, once its value is fixed we can detail the predictions of these 
energy loss models for the LHC: 

(i) Centrality dependence. 

Parton energy loss models that implement final state effects only, predict the 
absence of nuclear suppression in the most peripheral collisions. As a consequence, 
one expects that the centrality dependence of -RpbPb at the LHC parallels the one 
observed at RHIC. 

(ii) pj'-dependence of -RpbPb at the LHC. 

At RHIC energies, the slope of the partonic pT-spectrum gradually steepens as one 
moves from pt ~ 10 GeV/c to the absolute kinematic boundary pt = 100 GeV/c. 
This implies that to obtain the same value of Raa at higher px, a smaller fraction of 
parton energy loss is needed ("trigger bias effect", see reference [96]). In contrast, at 
the LHC, the partonic p-r-spectrum will show almost the same power-law over the 
entire range 10 GeV/c < px < 100 GeV/c, since this pr-range is far away from the 
kinematic boundary at the LHC. This implies that for a p^- independent Raa{pt)i 
one requires a constant, pr-independent fractional energy loss, not predicted in 
current energy loss models [99]. For a mild increase of e.g. gLHC — l-25gRHiC) 
models thus indicate that for central collisions, i?pbPb(10 GeV/c < pr < 20 GeV/c) 
is the same or slightly larger (by up to ~ 0.1) than at RHIC, and that -RpbPb(PT) 
increases gently by 0.1 — 0.2 from px = 10 GeV/c to pt ~ 100 GeV/iJ^. 

tt In this kinematic range, the partonic spectrum at the LHC is ghion dominated, while it is strongly 
quark dominated at RHIC. This difference in composition decreases the nuclear modification factor at 
LHC, compared to the one at RHIC, and compensates partially, but not completely for the trigger bias 
effect. 
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(iii) Dependence of i?pbPb on hadron identity. 

Light-flavoured hadron spectra are expected to show the same nuclear suppression 
independent of hadron species at sufficiently high pt > p^'^- If the particle species 
dependence at intermediate Pt < p^'^ is due to a medium-dependent hadronization 
mechanism (such as proposed for instance in recombination models [106-110]), then 
Pr'^i^/^NN) is expected to increase with -^/s^ by up to 2 — 3 GeV/c from RHIC 
to the LHC [110]. In contrast, if the scale p^"^ is mainly set by the time dilation 
of the hadronization time [111], implying that for pr > p^'^ hadronization occurs 
outside the medium, then one may expect that p^"^ is y/s^^ -independent . Thus, 
the ^s^jv -dependence of p^'^ may provide complementary information about the 
mechanism underlying the anomalous baryon-to- meson ratio at intermediate pt- 
For heavy-flavoured hadrons, parton energy loss models predict a hierarchy in the 
nuclear suppression [112-115], which can be characterized by "hcavy-to-light" ratios 
of the corresponding nuclear modification factors of heavy-fiavoured over light- 
flavoured hadrons [116]. For D-mesons, the charm mass is expected to be too small 
to contribute to a mass-dependent suppression above pt > 10 GeV/c. However, 
since light-flavoured hadrons at the LHC are dominated by gluon parents, the 
heavy-to-light ratio of D-mesons is sensitive to the colour charge dependence of 
parton energy loss and is expected to exceed unity by up to a factor ~ 1.5 in 
an extended p^-range. In contrast, for i?-mcsons, the mass effect is expected to 
put strong limits on medium-induced energy loss in an extended px-Tange. The 
resulting heavy- to- light ratios are expected to be significantly larger than for D- 
mesons, reaching a factor ~ 2 — 4, even for relatively low estimates of the value of 
q [116]. 

(iv) Single inclusive photon spectra. 

Based on a naive extrapolation from RHIC data, one may expect that photon 
spectra deviate only mildly from ~ 1 at the LHC. However, current parton 
energy loss models allow for mechanisms which may lead to significant medium- 
modifications: The medium-induced photon bremsstrahlung of hard partons may 
enhance the photon yield at high pt [117]; hard partons, which fragment into 
photons, may reduce the photon yield at high Pt [118]. Both effects are of order 
ftemj and it is unknown to what extent they cancel each other and how they may 
vary as a function of pj. 

7.1.2. Testing effects not encoded in current parton energy loss models. We now turn 
to phenomena that are qualitatively novel in the sense that they are not encoded in 
current parton energy loss models tested at RHIC, but may become important at the 
LHC. 

• Is high-pT hadron suppression at mid-rapidity a final state effect at all ^/s^l 
At RHIC, the absence of suppression in -RdAu showed conclusively, that initial state 
effects are unimportant for Raa [72-75]. However, high-pr suppression in pA has 
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Figure 6. The nuclear modification factor Raa as a function of transverse momentum 
at mid-rapidity. Data are for hadronic spectra measured at RHIC in ^s„„ = 200 GeV 
Au-Au collisions. The dashed line is a straight continuation of the high-p^ trend 
at RHIC, reproduced by models of radiative parton energy loss. Arrows indicate 
qualitative tendencies of how Raa may change at the LHC, see text for detailed 
discussion. 

been predicted for sufficiently high ^/s^ as a consequence of non-hnear QCD 
evolution in the so-called geometric scaling window [119-123]. At present, the 
value of ^/s^ at which non-linear evolution starts to become relevant for particle 
production at intermediate pt is unclear. If the onset of such non-linear evolution 
effects should lie below LHC energy, then this would be signalled by i^ppb < 1 at 
LHC mid-rapidity. The enhancement i^dAu > 1 observed at RHIC mid-rapidity 
would turn into a characteristic suppression as a function of ^/s^. Moreover, one 
expects in this case that the nuclear modification factor in Pb-Pb will be suppressed 
due to both, initial and final state effects. 

• Does RphPhiPr) show indications of medium- dependent Q^-evolution? 

QCD evolution underlies the partonic fragmentation process in the perturbative 
regime, as well as the scale dependence of fragmentation functions. The question of 
how the medium affects the QCD scale evolution is difficult to address theoretically, 
since medium effects are "higher twist" , i.e. subleading by powers of Q^, 
although they may be nuclear enhanced by geometric factors oc A^^^ [124, 125]. 
Qualitatively, however, one may expect that — despite the time dilation of 
the parton fragmentation process in the target rest frame — at sufficiently 
large Q^, parton splitting occurs on length scales too short to be resolved by 
the medium, and short- distance contributions to parton fragmentation remain 
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unmodified accordingly. It has been speculated that this may lead to a significant 
rise of Raa{pt) in the range 20 GeV/c < pr < 100 GcV/c [126], but detailed 
model studies are still missing. More generally, the fact that LHC has access to a 
logarithmically wide p^-range may provide novel opportunities to test the medium 
dependence of QCD-evolution. 

• What is the dynamical mechanism underlying the nuclear modification of 
quarkonium? 

Data on R'^^{pt,7]) at RHIC are known up to less than px < 5 GeV/c. They 
indicate suppression at small pt, Raa{Pt ^ IGeV /c,\r]\ < 0.35) ~ 0.3 — 0.4, 
and possibly an increase of R"^^ with pt- This is qualitatively similar to i?^^ 
of light hadrons. However, the physics invoked to account for the suppression 
of J/^'s shows marked differences if compared to the high-p^ suppression of 
light hadrons. Quarkonium suppression arises from the fact that the attraction 
between heavy quarks and anti-quarks weakens with increasing temperature due 
to dynamic screening effects [127]. According to recent models, directly produced 
J/\l''s may not dissociate until well above the energy densities attained at RHIC, 
but the Xc's and ^'"s, whose decay contributions are estimated to ~ 60% of the 
J/^ yield in hadronic coUisions, are expected to dissociate. The observed value 
Raa(pt < 1 GeV/c, \r)\ < 0.35) ~ 0.3 - 0.4 thus appears natural, if one assumes 
that all excited cc bound states are dissociated in the medium [128, 129]. 
At the LHC, one will have for the first time experimental access to significant rates 
of both bound charmonium and bottomium states. The lowest lying 66-states are 
more tightly bound than the charmonium states, and one thus expects that they 
dissociate at higher temperature. Accordingly, in dissociation models one expects 
quite generally that the nuclear suppression of these bottomium yields is smaller or 
at most as large as that of the corresponding charmonium yields. A qualitatively 
opposite behaviour can be expected, if secondary production mechanisms, such 
as recombination, start playing a role [23]. These become more effective with 
the number of heavy quarks produced, and are thus more efficient in enhancing 
charmonium bound states. 

At high transverse momentum, pt>5 — 7 GeV/c, recombination effects arc absent, 
and LHC will establish how the high-pr modification of the spectra of mesonic 
QQ bound states differs from that of heavy-light and light-light fiavoured mesons. 
The following quahtative considerations illustrate the importance of formation time 
effects in this context: i) If high-p^ J/^'s originate from gluon parents and if 
these gluons propagate over long distances before hadronizing, then the nuclear 
modification factor of J/^^'s is expected to be reduced due to gluon energy loss, 
and should take the same value as that for light-fiavourcd hadrons. ii) In contrast, 
if high-po" J/\E''s would originate from the fragmentation of c-quark parents and 
if the quark propagates over long distances prior to hadronizing, then the nuclear 
modification factor is expected to be reduced due to quark energy loss only, and 
should match the smaller reduction of heavy-light fiavoured hadrons R"^^ ~ i?^^ > 
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R\j^- In this case, one would also expect in pp collisions open charm production 
associated with J/\E' production, iii) If in contrast the formation time of the J/\E' 
is small compared to the extension of the collision region ~ 10 fm, then it is a QQ 
bound state of velocity v increasing with p^, which propagates through the medium. 
For this case, non-perturbative calculations based on the AdS/CFT correspondence 
suggest that the screening length Lg of the QQ potential at temperature T shows 
a characteristic velocity scaling Lg{v,T) oc Ls{0,T)/^y^ [130]. Model estimates 
indicate that depending on the binding energy of the QQ-state, it is this velocity- 
dependent dissociation effect which may dominate in an extended intermediate 
PT-regime. Only in this latter case will different bound states such as J/ ^ and ^' 
show a different degree of suppression. 

• Are high-pT single inclusive photon spectra sensitive to final state medium effects? 
In hadronic collisions at the LHC, a significant fraction of the single inclusive pho- 
tons arises from the fragmentation of quarks and (at next-to-leading order) gluons. 
If the parent partons suffer medium-induced energy degradation prior to fragment- 
ing into a photon [118], then i?^^ is reduced and the elliptic flow of photons receives 
a positive contribution [131, 132]. The strength of this effect may allow one to con- 
strain the photon formation time but is difficult to estimate a priori. On the other 
hand, the interaction of the produced partons with the medium can lead to addi- 
tional bremsstrahlung photons [117]. This effect increases -R^^, and contributes to 
a negative V2 for photons [132]. To the best of our knowledge, these expectations 
remain qualitative to date, but the wider kinematic reach at LHC should help to 
disentangle them in the data. We flnally note that similar modifications of the 
high-pT spectrum are not expected for Z-bosons, since all aspects of its production 
are local and unlikely to interfere with the typical length and momentum scales 
present in the produced QCD matter. 

7.2. The nuclear modification factor as a function of rapidity. 

7.2.1. Evolution of the Cronin peak at RHIC. The notion "Cronin peak" refers to 
the enhancement of the nuclear modification factor R^a above unity at intermediate 
transverse moment, 2 GeV/c ^ Pr ^ 4: GeV/c. At lower transverse momentum, the 
nuclear modification factor is generally suppressed RhA < 1- We hasten to remark that 
the notion "Cronin peak" is but a description and not an explanation of the shape 
of RhA- As discussed below, there is no completely satisfactory explanation of the 
dynamical origin of the Cronin peak so far, but LHC is well-positioned to provide 
additional insight. 

The kinematically available p^-range decreases with increasing rapidity, and data 
from RHIC are currently hmited to pt < 5 GeV/c for I77I > 2. Within this hmited 
range, one observes the following apparently generic trends in RHIC data [133] : 

(i) Rapidity and pr-dependence of RdAu- 
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At RHIC, -RdAu shows a typical Cronin enhancement for rj < 1. This Cronin peak 
in the range 2 GeV/c < Pt ^ 4 GcV/c monotonically decreases in the deuteron 
fragmentation region. One finds i?dAu — 1 at ~ 1, -RdAu — 0.5 at 77 = 3 [134-136] 
and even smaller values for 7r°'s at higher rapidity [137]. In contrast, in the 
nucleus fragmentation region, the ratio Rep of central over an equivalent number of 
peripheral collisions, which is closely related to RdAu, shows an enhancement [136] 
above unity, which increases with increasing rapidity. 

(ii) Centrality dependence of suppression pattern in d-Au. 

At mid-rapidity, it!cp increases around the Cronin peak for increasing centrality. In 
contrast, the centrality dependence at forward deuteron rapidity is inverted [136]: 
-RdAu decreases with increasing centrality for 77 > 1, and the centrality dependence 
at 77 ~ 1 is neghgible. 

In the context of RHIC data, the discussion of these phenomena has focused on two 
initial state effects, which we address now: multiple scattering and non-linear QCD 
evolution of the incoming parton distribution functions. 

The Cronin peak is often thought of as the consequence of a multiple scattering 
picture, in which the partons in the deuteron wavefunction undergo multiple interactions 
in the target nucleus prior to producing relatively high-py hadrons. Incoherent multiple 
scattering of these incoming partons leads to a transverse momentum broadening of 
the initial parton distribution, which translates into a correspondingly broadened single 
inclusive hadron spectrum. This can account for the observed Cronin peak at mid- 
rapidity [138, 139], though it is unclear whether it can account for the particle species 
dependence of the effect. Moreover, at least in their current model implementations, 
such multiple scattering models predict the persistence of the Cronin peak at forward 
rapidity, and they imply an increase of p^-broadening with increasing centrality at all 
rapidities [138]. This contradicts the generic trends seen in the RHIC data. 

In studies of non-linear QCD evolution at small momentum fractions x, one 
gcncrically finds that the growth of unintegrated gluon distribution functions with 
In 1/a; (or with In ^/s^) is saturated up to a scale pt < Qs {x) , which grows with 
1/x. Moreover, above this saturation scale, non-linear QCD evolution characteristically 
changes the power-law in a wide geometric scaling window. Contact between 
these findings and the phenomenology of d-Au collisions is made by the following 
observations: First, toward deuteron projectile rapidity, smaller momentum fractions 
Xau of the nuclear parton distribution functions become relevant for particle production. 
Eventually, xau will become small enough for non-linear QCD evolution to be applicable. 
Second, convoluting non-linear evolved unintegrated nuclear parton distributions 
schematically with hard processes, one finds that non-linear QCD evolution implies 
the decrease of RhA with increasing Inl/xAu [119-123]. This provides a conceivable 
explanation for the rapidity dependence of RdAu- 

Non-hnear QCD evolution does not account for all trends seen in the data. In 
particular, the Cronin peak itself is not a dynamical consequence of non-linear QCD 
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evolution, it is just a conceivable initial condition, which is quickly washed out by the 
evolution. So, if saturation physics is the correct explanation for the rapidity dependence 
of -RdAu, then one knows that it is not applicable at RHIC mid-rapidity. Also, while 
particle species identified data on -RdAu are not accurate enough to allow for decisive 
tests, the question whether a purely partonic explanation is sufficient to account for 
i?dAu in the experimentally tested range remained open so far. Moreover, it is unclear 
on theoretical grounds whether saturation physics can be expected to apply for the 
relatively large values xau > 0.02 [140] which dominate forward particle production at 
RHIC. 

Experiments at the LHC will allow us to compare the ^/s^- and 77- dependence. 
To illustrate that this may provide a decisive test for current models, let us consider 
an alternative explanation of the rapidity dependence of -RdAu at RHIC, based on the 
following picture: Partons in the deuteron wave function undergo inelastic multiple 
scatterings on the nuclear target field. Hence, these partons split due to interactions 
with the target. Since splitting is very effective in energy degradation, this will deplete 
the hadron yield at forward rapidity but will enhance it at mid-rapidity. Also, the 
opposite centrality dependence at mid-rapidity and forward rapidity [point (ii)] can be 
understood in this way. In this picture, the suppression in -RdAu increases strongly 
with increasing projectile rapidity, because the energy degradation occurs on top of 
an increasingly steeply falling spectrum at forward rapidity. Transverse momentum 
broadening could still contribute to the Cronin peak at mid-rapidity, but would not be 
able to overcome the reduction at forward rapidity. This picture accounts for the same 
rapidity dependence as saturation models. In contrast to saturation models, however, it 
implies that the Cronin peak will not disappear with increasing ^/s^, but will persist 
at LHC energies. 

7.2.2. Conceivable effects on the t]- dependence of RpPh cind -RpbPb o-t the LHC. The 
example given above illustrates that comparing the y/s^- and //-dependence at the 
LHC will provide a qualitatively novel test for the saturation physics interpretations of 
measurements at RHIC. This is so, since an increase in both y/s^^ or t] gives access 
to smaller momentum fractions x in the parton distributions, and thus has similar 
implications in saturation models. At the LHC, the ^s^^^^ -dependence of measurements, 
in combination with their //-dependence, will become a tool to discriminate effects from 
small-x QCD evolution from other conceivable mechanisms. 

Focusing in the following solely on the nuclear modification factor, we now list 
phenomena which may affect significantly the r^-dependence of -RpbPb at the LHC, and 
which may be disentangled by studying the ^s^^^ -dependence in a combination of data 
from RHIC and LHC. Our list starts with conceivable final-state effects: 

(i) Dependence of -RpbPb on dNch/drj for fixed centrality. 

As discussed in section 17.1.1^ the quenching parameter q is expected to grow 
monotonously with dNch/drj. So, the quenching parameter should be smaller at 
forward rapidity. This effect contributes to an increase of -RpbPb with rj. 
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(ii) The r/- dependence of the partonic pr-spectrum. 

With increasing 77, partonic p^-spectra get steeper. This is a simple kinematic 
effect, present both at RHIC and at the LHC, but quantitatively different. As a 
consequence of this effect, the same amount of parton energy loss leads to a decrease 
of -RpbPb with rj (trigger bias). 

(iii) Flow effects on i?pbPb- 

The initial parton, produced in a hard collision, needs not be produced within 
longitudinal comoving matter. In case that it is not, there is a relative longitudinal 
velocity between the hard projectile and the medium, and energy loss is expected to 
be higher. This effect is hkely to contributes to a decrease of -RpbPb with 77 [141,142], 
though estimates of its magnitude vary widely [143] . 

(iv) Initial state effects. 

As discussed above, the nuclear modification of parton distribution functions is 
expected to affect -RpbPb- In particular, models based on non-linear small- a; 
evolution predict [119-123] an additional decrease of i?pbPb with 77. These models 
can be tested by comparing the ^/s^- and ?]- dependences in p-Pb collisions. 

8. "Jet-like" particle correlations and jets 

The leading hadronic fragments of highly energetic parent partons, measured in single 
inclusive hadron spectra, are strongly modified at RHIC and they are expected to be 
strongly modified at the LHC. Any model of the dynamical mechanism underlying this 
medium modification has implications for the entire parton fragmentation pattern, and 
that is: jets and jet-like observables. Jet measurements in heavy ion collisions are 
sensitive to how high-energy partons are attenuated in matter and how they equilibrate 
kinetically and chemically. In turn, these medium-modifications of jet fragmentation 
characterize properties of the produced medium. 

In general, parton fragmentation leads to multiplicity distributions with broad 
variances. As a consequence, any particle trigger used to select jet-like observables will 
bias significantly the fragmentation pattern. Even prior to invoking medium effects, 
such biases have dramatic consequences: In a typical single inclusive hadron spectrum 
(i.e. single particle trigger) at px > 20 GeV/c, the hadrons will typically carry on 
average ~ 3/4 of the energy of their parent partons. In contrast, the leading hadron in 
a 100 GeV/c jet, initiated by a light parton, carries typically only ~ 1/4 of the jet energy, 
simply because this jet fragmentation pattern is not biased by a single particle trigger. 
In the presence of a medium, additional "trigger biases" may arise. For instance, in the 
presence of strong final state energy loss, a high-p^- particle trigger will select particles 
produced mainly at the outskirts of the nuclear overlap region. The parent partons of 
these hadrons have had a particularly small in-medium path length and thus suffered 
particularly little parton energy loss (surface bias) [144]. Also, a high-p-r particle trigger 
will prefer events in which the initial state p-r-broadening effects move the dijet invariant 
mass towards the trigger. Thirdly, triggering on a high-energy hadron or requiring 
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a jet can lead to structures in the distribution of soft "background" particles, which 
are typically counted towards the medium, but which are related to the trigger and 
would not be found in minimum bias events. These general considerations prompt us to 
distinguish in the following discussion between "true" jets, jet-like particle correlations 
and soft structures causally related to high-pr triggers. 

8.1. The medium-modification of "true" jets 

"True" jet measurements, that is measurements of the entire fragmentation pattern of 
high-£^T parent partons, have not been performed in heavy ion collisions so far. In 
the context of RHIC data, "jet quenching" refers to the suppression of single inclusive 
hadron spectra and high-p^ particle correlations. Yet, measurements at RHIC, as well 
as models of parton energy loss, give rise to a set of general expectations for "true" jet 
measurements in heavy ion collisions: 

(i) Longitudinal jet multiplicity distributions soften. 

Parton energy loss, combined with energy-momentum conservation implies that 
the energy lost by the leading parton or hadron in the parton shower reappears in 
additional multiplicity of softer fragments. The entire longitudinal jet multiplicity 
distribution is expected to soften, and the total jet multiplicity is expected to 
increase, sec e.g. references [126,145]. 

(ii) Transverse jet multiplicity distributions broaden. 

Essentially all models of parton energy loss assume a significant transverse 
momentum transfer from the medium to the jet projectile. As a consequence, 
parton energy loss is generally thought to be accompanied by a broadening of the 
jet fragmentation pattern in the plane orthogonal to the jet axis [141, 142, 146, 147]. 
In case that the momentum transfer from the medium is asymmetric, for instance 
since the parton is embedded in a collective flow field, this jet broadening may show 
characteristic asymmetries [142,146]. 

(iii) The hadrochemical composition of jet fragments may be modified. 

To date, most studies of jet medium-modifications focus on the longitudinal and 
transverse energy and multiplicity distributions. However, in current models of 
parton energy loss, the medium couples to the parton shower via gluon exchange, 
and thus alters the colour flow in the shower. This may be expected to affect the 
hadrochemical composition of the jet. Also, in principle, other quantum numbers 
such as baryon number or flavour may be exchanged between the medium and the 
jet [148]. 

We note that even if the average longitudinal jet multiplicity distribution softens, 
it may be possible that high-p^ triggered particle correlations remain insensitive to the 
properties of the medium. This is so, for instance, if the high-pr trigger should select 
the subset of parton fragmentation patterns, that escaped with a negligible medium 
modification e.g. due to a surface bias effect. Similar remarks apply to the transverse 
jet multiplicity distribution and hadrochemical composition. 
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8.2. Jet-like particle correlations and a potential all-or-nothing mechanism 

There is a class of measurements, in which a trigger hadron of high transverse momentum 
p^'^ is correlated with associated hadrons as a function of their transverse momentum 
p^°'^ and their azimuthal angle A0 with respect to the trigger particle. We call such 
correlations "jet-like", ii p?p^°^ is relatively large, 2 GeV/c < p^^"'^ < p^'^, say. For a 
first theoretical work on jet-like correlations, see e.g. [149]. The generic trends seen in 
such correlation functions at RHIC are: 

(i) Near-side jet-like particle correlations in Au-Au are independent of centrality and 
similar to those in pp or d-Au. 

In pp and Au-Au collisions at RHIC, near-side (i.e. small A0) two-particle jet- 
like correlations show an enhancement characteristic of hard-scattering processes. 
Compared to pp collisions, the yield of high-pT trigger particles decreases by a 
factor ~ 5 from peripheral to central Au-Au collisions at RHIC. In contrast, jet- 
like two-particle correlations do not show a significant centrality dependence. For 
sufficiently high threshold trigger 8 GeV/c < p^'^ < 15 GeV/c, the yield and 
A(/)-width of the near-side distribution is insensitive to the centrality of Au-Au 
collisions, and coincides with the measurement in d-Au collisions [150]. The same 
has been observed for lower trigger thresholds [151]. Also other features of jet- 
like pT-triggered correlation functions, such as the ratio of like-sign to unlike-sign 
pairs in jet-like correlations [151], do not show any centrality dependence and are 
consistent with the data found in pp collisions. 

(ii) Back-side jet-like particle correlations decrease in yield with increasing centrality, 
but keep approximately the same width. 

For intermediate pt triggers (4 GeV/c < p^'^ < 6 GeV/c) at RHIC, the associated 
particle yield for p^^°^ > 2 GeV/ c disappears as a function of centrality [151] . If one 
raises the trigger threshold to higher values (8 GeV/c < p^'^ < 15 GeV/c), then the 
back-side jet-like structure reappears again, but the yield strongly decreases with 
centrality. The back-side structure shows no sign of azimuthal broadening [150, 152] . 

The above features are qualitatively consistent with a schematic all-or-nothing 
mechanism, based on the following picture: If a hadron is triggered on with a high 
p^'^, then it is the leading fragment of a parton shower, which propagated essentially 
unperturbed through the medium ("complete survival of entire jet structure"). On the 
other hand, if the parton shower is significantly perturbed by the medium, then the 
energy of the leading fragment is degraded to such an extent, that it becomes unlikely 
to find this fragment in a high-p^- trigger bin ( "no survival at all" ) . This all-or-nothing 
picture may be regarded as the most extreme form of a trigger bias: the trigger selects 
the subclass of unmodified parton fragmentation patterns and the medium-modification 
establishes itself solely in the reduced yield. In this way, this all-or-nothing picture 
accounts for the strongly reduced yield of high-pT trigger particles, characterized e.g. 
by the nuclear modification factor, as well as for the suppression of the back-side yield. 
It can also account for the absence of broadening in both the near-side and the away- 
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side peaks by arguing that the particle pairs entering the jet-hke correlation function 
belong to parton showers which escaped the medium essentially without interaction and 
thus without signs of medium-induced broadening. The picture is also qualitatively 
consistent with finer features seen in the data, such as the observation that for near- 
side correlations, the particle yield as a function of the effective fragmentation variable 
zt = p^'^^'/p^^'^ is the same in d-Au and Au-Au, independent of centrality; on the away- 
side, the particle yield decreases with centrality but shows the same z-p-slope [152]. 

Can such an all-or-nothing mechanism be consistent with the dynamics of QCD 
radiation physics? To address this question, one may note first that for a steeply falling 
partonic pr-distribution, it is conceivable that all high-pr trigger bins are dominated 
by hadrons, whose parent partons suffered no medium- induced parton energy loss [96]. 
In other words: while hadrons, whose parents suffered some medium-induced energy 
loss must end up in some p^-bin, they can — for steeply falling distributions — 
always be shifted to an abundantly populated lower p^^-bin, in which their yield is 
statistically negligible. Recent implementations of radiative parton energy loss can 
account at least qualitatively for this possibility by two features [98,99,102]: first, 
even for dense systems, recent models allow for a sizeable finite probability that the 
parton shower propagates unperturbed through the medium. Second, the distribution 
of leading fragments in the parton shower turns out to be very fragile, once the parton 
shower has interacted with the medium. In this way, current model implementations 
contain the main ingredients needed for implementing a strong surface bias, which may 
underlie the all-or-nothing mechanism sketched above. In model studies, one has also 
addressed more refined questions, such as whether the surviving yield in the away-side 
correlation arises predominantly from particle pairs emitted tangentially to the surface 
of the collision region [153,154], so that neither the trigger nor the associated recoil 
particle traverses a significant amount of matter. 

The all-or-nothing mechanism outlined here is a working hypothesis, which finds 
some support in RHIC data and current model analyses. If true, it is a dramatic 
illustration that jet-like particle correlations fall short of characterizing the distributions 
of quenched jets, simply because they trigger mainly on the small fraction of unquenched 
survivor jets. To refine this all-or-nothing mechanism (or rather: to replace it by 
a picture which allows for gradual manifestations of parton energy loss on jet-like 
correlations), one should study in particular correlations with lower p^'^°'^. This is 
so, since the trigger particle of p^'^, to the extent to which it does not arise from a 
medium-independent fragmentation, should be accompanied by an increased associated 
yield at sufficiently small p^^°^. At RHIC, lowering p^^^oc i^gjQ^ 2 GeV/c for Au-Au 
collisions, one has observed indeed an enhanced associated yield with clear indications 
of broadening of the away-side peak. However, the kinematic range p^^"'^ < 2 GeV/c 
is difficult to disentangle from the large underlying event multiplicity and it may be 
affected by other mechanisms, see section 18.31 below. 

A trigger-biased class of jet measurements, which shows medium-modifications 
of associated jet multiplicity only below p^^°^ < 2 GeV/c, provides arguably only 
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limited access to a study of the entire quenched jet fragmentation — except, of course, 
if one could demonstrate that this trigger bias is unimportant and that these jet- 
like correlations are characteristic for the average medium- modified parton shower. 
The wider kinematic reach of heavy ion collisions at the LHC may provide means 
to this end. A jet of = 200 GeV has on average ~ 7 charged hadrons with 
^assoc y 5 GeV/c. Although jet-like correlations based on single trigger particles will 
bias significantly the average jet fragmentation pattern, one expects qualitatively that 
the distribution of associated particles should show imprints of medium-modifications 
(namely signs of pr-broadening and enhanced yield) in a wider range of p^^°^, which 
can be disentangled more clearly from the underlying event multiplicity. However, this 
qualitative expectation is not yet supported by model studies. 

8.3. The pedestal, the ridge, the Mach cone and all that ... 

High-pT triggers affect the underlying event in hadronic collisions. For instance, in 
comparison to minimum bias data, triggering on a high-pj- hadron in a pp collision 
increases the soft event multiplicity by a factor of order ~ 2. The hard parton sits 
on top of a "pedestal", which is wide in rapidity [155]. Within perturbation theory, 
such a pedestal may be expected, since large Q^-processes are accompanied by initial 
state radiation, which is broad in rapidity and which will manifest itself in additional 
low-pT hadrons. Since this initial state radiation moves over long ranges with the beam 
fragments, non-perturbative physics may play an important role as well. Also, multiple 
parton interactions may contribute to the pedestal effect [156] . 

The pedestal observed in high-pr-triggered hadron collisions is the prototype of a 
phenomenon, which is clearly related to the presence of a high-Q^ process, but which 
is not due to final state parton fragmentation. As such, the pedestal is a structure in 
the low-pT trigger-associated particle yield, which one cannot expect to reproduce in a 
model that superimposes a high-p^^ final state fragmentation pattern on the multiplicity 
distribution of a minimum bias event. The state of the art of modeUing high-p^ 
phenomena in heavy ion collisions is of the latter type, and one wonders whether there 
are — like the pedestal in pp — characteristic features in heavy ion collisions, which 
one misses in models superimposing medium-modified hard processes on minimum bias 
soft background. 

One candidate for such a feature is the "ridge" : a trigger particle is accompanied 
by additional associated hadronic activity in some range of intermediate p^^°'^ at the 
near-side only. This additional multiplicity is wide in rapidity but, unlike the pedestal, 
it is not balanced by a similar amount of activity in the same range of p^^°^ on the 
away-side. This phenomenon may arise e.g. in a picture [157], in which the pedestal 
is embedded in a transverse fiow field. Namely, triggering on a high-pr particle, one 
selects an interaction point which will preferably lie away from the centre of the collision 
region towards the direction of the trigger pr. At this point in the transverse plane, 
collective transverse fiow is also expected to point in the direction of the trigger p^. 
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So, any additional initial state hadronic activity, associated with this trigger, may be 
expected to be transported by transverse flow towards the near-side. 

The above is but one, albeit speculative, illustration that if one aims at studying 
values of p^^°'^ comparable to those in the bulk multiplicity, the study of medium- 
modified jet measurements cannot be limited to the study of medium-modified final state 
parton fragmentation patterns on top of minimum bias events. For lowp^^^"*^, it becomes 
difficult to establish which part of the additional hadronic activity emerges from the 
fragmentation and energy loss of a hard final state parton. We note that also the much- 
discussed, broad structures in the away-side correlations, which have been suggested to 
indicate the appearance of Mach cones, do not persist for higher but are only 

seen in a rather narrow range of low transverse momentum. Radial fiow, anisotropic 
flow, initial state radiation and trigger bias effects may all affect characteristic features 
of associated particle distributions in this low p^^°'^ -regime. Heavy ion collisions at the 
LHC may help to clarify the dynamical understanding of such soft structures related 
to high-pT trigger particles, since the hadronic activity in both the incoming and the 
outgoing state is expected to increase significantly with the trigger p^^, and may manifest 
itself in a wider range of p^°'^. 

9. Connecting Heavy Ion Phenomenology with QCD 

How collective phenomena emerge from the fundamental laws of elementary particle 
physics is a multi-faceted question, which in the range of extreme matter densities, 
where physics is determined by partonic degrees of freedom, has been addressed 
in different theoretical approaches. Historically, high-temperature QCD equilibrium 
dynamics, studied non-perturbatively in lattice calculations or perturbatively in finite 
temperature field theory, has been the first theoretical approach with the potential of 
connecting heavy ion phenomenology with first principles of QCD. In particular, the 
most dramatic collective phenomenon, expected in finite temperature QCD, namely the 
phase transition to a quark gluon plasma at a critical temperature and baryochemical 
potential, has been firmly established in lattice QCD. By now, these techniques 
are applied to many questions of phenomenological relevance at the LHC, including 
quarkonium suppression, the medium-modification of spectral functions, dissipative 
transport coefficients, and fiuctuation measurements [158,159]. 

On the other hand, heavy ion phenomenology has estabhshed over the last 
two decades strong indications that effects of directed collective motion are at least 
as important for understanding the dynamics of heavy ion collisions, than effects 
of random thermal motion. These two concepts, collective dynamics and local 
equilibrium, can coexist. Indeed, the modelling of heavy ion collisions in terms 
of perfect fluid dynamics illustrates the extent to which a mesoscopic system with 
extreme position-momentum gradients may still maintain local thermal equilibrium. 
In the discussion of this hydrodynamic modelling, the emphasis has shifted gradually 
from fundamental tests of QCD thermodynamics (namely whether the QCD phase 
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transition and its order leave traces in the dynamical evolution) to fundamental tests 
of QCD hydrodynamics (namely the test of dissipative properties of the matter, such 
as viscosities). There are many reasons for this gradual shift of focus at collider 
energies, starting with the observation that at RHIC and LHC energies one likely 
overshoots the critical energy density significantly, and ending with the notorious 
problem of identifying "unambiguous" signatures of the QCD phase transition. In 
principle, the characterization of hydrodynamic features in heavy ion collisions provides 
an opportunity of connecting heavy ion phenomenology to first principles in QCD, since 
both properties of the QCD phase transition and dissipative transport coefficients are 
calculable directly from the QCD Lagrangian. In practice, however, one prerequisite 
for exploiting these opportunities is a very good experimental and theoretical control 
over the "perfect fluid baseline" on top of which one aims at establishing dissipative 
properties. Our discussion in sections [3l [5] and [6] also identifled how measurements at 
the LHC can help to establish whether conditions close to this perfect fluid baseline are 
realized at the LHC and whether they were realized at RHIC. 

It would be an unwanted bias to limit the study of "soft" physics at the LHC to 
manifestations of QCD thermo- and hydrodynamics. In this review, we deliberately 
started from the observation that several apparently generic trends in the existing data 
(e.g. in multiplicity distributions and collective flow) have not yet found a satisfactory 
explanation. Agnostic extrapolations of these trends to the LHC are at odds with 
the extrapolation of current models, be it hydrodynamics or saturation physics. This 
indicates that LHC will be a discovery machine also in the soft physics sector. In 
particular, the support for an interpretation of data in terms of hydrodynamics or 
saturation physics would be strengthened qualitatively, if one discovered at the LHC 
deviations from the so far apparently generic trends, which are characteristic for the 
currently advocated dynamical models (such as a mild but distinct power-law y/s^ 
increase of event multiplicity, or a deviation of V2 from In ^s^-scaling) . On the other 
hand, a conflrmation of these trends may prompt us to reassess our understanding of 
the soft matter produced in heavy ion collisions at both LHC and RHIC. 

LHC will also be a discovery machine outside the soft physics sector at mid-rapidity. 
This is so mainly because of the logarithmically wide range in transverse and longitudinal 
momentum, which opens up at 30 times higher ^/s^ = 5.5 TeV. As discussed in 
sections [7] and [H the resulting abundance of hard processes at the LHC provides many 
novel tools for probing the produced soft matter. The prerequisite for exploiting this 
opportunity is a very good experimental and theoretical control over how the medium 
modifies hard processes due to interactions in the final and initial state. Our discussion in 
section[8]identified how measurements at the LHC can improve this control, in particular 
by extending jet quenching studies significantly beyond the analysis of medium-modified 
leading fragments. This is likely to refine our understanding of hard probes at the LHC 
and at RHIC. 

It would be an unwanted bias to limit the study of hard physics at the LHC to its use 
as "hard probes" . With the significantly wider kinematic reach, heavy ion physics gains 
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experimental access to other fundamental properties of QCD. In particular, medium 
modifications of the QCD-evolution in both (mainly via the transverse momentum 
dependence) and In l/x (mainly via the rapidity and ^s^^^^ -dependence) become testable 
at the LHC. Again, characteristic deviations of LHC measurements from the agnostic 
extrapolations discussed here may provide some of the cleanest possibilities of identifying 
and ultimately quantifying the manifestations of medium-dependent QCD evolution. 
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